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Abstract. Many recent experiments employ several parametric down conversion 
(PDC) sources to get multiphoton interference. Such interference has applications 
in quantum information. We study here how effects due to photon statistics, 
misalignment, and partial distinguishability of the PDC pairs originating from different 
sources may lower the interference contrast in the multiphoton experiments. 

PACS numbers: 03.65.Ud, 03.67.Bg 
1. Introduction 

Optical setups are an important implementation of various quantum information tasks, 
e.g. quantum teleportation P, entanglement swapping [2] or quantum repeaters [3]. 
Most of the modern quantum optical experiments involving few photon interference use 
spontaneous parametric down-conversion (PDC) sources. It is a highly unlikely process 
in a non-linear crystal, in which a photon from the pumping beam changes into two 
entangled photons in the output modes. The quantum correlations are manifested in 
polarizations, momenta, and frequencies of the emitted pairs '5]. 

Recent years have brought a significant progress in realizing multisource PDC 
experiments allowing multiphoton interference. It was possible to experimentally verify 
five-photon entanglement 0, and recently, six-photon interference was observed [7j. 
However, these results required a considerable effort, putting in question the feasibility 
of realizing more complex systems. Therefore it is important to check if there are 
any fundamental obstacles for conducting future multiphoton interference experiments. 
Such obstacles might make some involved quantum information protocols impossible to 
realize via purely photonic techniques. 
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Visibility (interferometric contrast) in multiphoton PDC experiments may be 
impaired by statistical properties of the emission process p], misalignment and partial 
distinguisliability of the entangled pairs. Therefore, a quantitative analysis of the 
influence of these factors is important, if we want to use non-classical properties of 
light in quantum communication. E.g., in order to show a conflict between quantum 
mechanics and classical description, the visibility should be greater than some critical 
value V'^"^. Different critical values might be needed to show that the state has 
entanglement of a required type. This report investigates what is the impact of the 
specific properties of PDC on the observed multiphoton visibility, and what limitations 
follow from the necessity to breach these values. 

We calculate the maximal possible visibilities under influence of these factors. 
Therefore, the description will be as simple as possible, as inclusion of other traits 
of the experiments would work toward lowering of the achievable visibility. To study 
the statistical properties we shall use the simplest description of type-II parametric 
down-conversion with just four modes of the radiation. We shall ignore completely the 
frequency-momentum structure of the states of the emitted signal-idler pairs. However, 
when studying effects due to distinguishabihty of such pairs originating from different 
emission acts, we shall ignore the statistical properties, and derive our results basing 
solely on the structure of the photon pair states. The other principal approximation that 
we shall make is neglecting in our description all sources of losses. Thus the problems 
of detection and collection efficiency will be ignored. 

The paper is organized in the following manner. In section II we study the decrease 
of visibility due to the production of the additional pairs of photons, which occurs in the 
strong pumping regime. This will be discussed for a Bell-type experiment and a setup 
designed to observe GHZ-correlations. Section III is divided in two parts. The first 
analyzes the problem of mode mismatch, which may occur if the interfering photons 
come from two different sources. In the second part we first give a simple description 
of the spectral properties of the effective two-photon output state of the spontaneous 
parametric down conversion process. This is followed by a consideration of the reduction 
of interferometric visibility due to the strong temporal correlations, which characterize 
a such a state. This affects the indistinguishability of the quantum processes associated 
with different emission acts, and therefore interference. The usual method of enforcing 
indistinguishability by employing suitable filters is studied quantitatively for a sequence 
of multiphoton experiments. We start with a single source experiment (two photon 
interference) and continue the analysis up to five sources (ten photon interference). In 
section IV we analyze the properties of the effective output state in PDC process. The 
basic result of this section is an observation that the reduction of the interferometric 
contrast due to (unwanted) multiple emissions, in some interesting cases, cannot be 
modeled by an additional white noise admixture to the observed interference pattern. 
The noise in some cases has a much more involved structure. We also give a new 
definition of an overall interferometric contrast, which seems to describe the situation 
much better than the standard one. It involves interferometric processes of various 
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Figure 1. The scheme for the polarization analysis of entangled pairs. PDC - 
parametric down-converting source, PA - polarization analyzer (the measurement 
settings are parametrized by 0^, a^ - PDC output modes, A^-i- - analyzer output 
channels, Di± - detectors. 



order. Section V contains conclusions and summary. 

2. The influence of statistics for high conversion efliciency 

To obtain multiphoton interference with several PDC sources one has to go to sufficiently 
high efficiencies of the down conversion process, which is usually described as a 
"strong pumping regime" , because otherwise losses and detection inefficiency make the 
multiphoton coincidences prohibitively low. However, multiple emissions cause effects, 
which lower the visibility of the interferometric processes. 

2.1. Statistical properties of the PDC radiation 

Assume that the PDC source emits pairs of polarization entangled photons into the 
spatial modes ai and a2 (see Fig. [1]). The most simplified form of the Hamiltonian for 
the PDC process can be put as 

^ = ^x(«I,//«2,/f + + h.c. (1) 

Horizontally (H) and vertically (V) polarized photons occupy two spatial modes (oi, 02), 
and X is an effective coupling constant proportional, among others, to the pump power 
and the second-order nonlinearity of the crystal. After the interaction time the resulting 
photon state is given hj = e~*^*|f2), namely, 

^ 00 

1^) = E tanh"+-ir|n)fjm)rjn)fjm)r„ (2) 

n,m=0 

where \Q) is the vacuum state, K = xt and \i)xy denotes i-photon state in the spatial 
mode X and polarization mode y, see e. g. ||9]. We shall now study the interference 
effects obtainable when one suitably overlaps radiation of several such sources. 

2.2. A Bell-type experiment 

Let us first discuss a standard Bell-type experiment involving just one source. One 
can pass the radiation given by eq. ([2]) through polarization analyzers PAi and PA2 
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Figure 2. The visibility of two photon interference in the experiment of Fig. 1 as the 
function of the PDC process efficiency parameter K. 



(Fig. 1). The polarization measurements are performed by operators corresponding to 
qubit dichotomic observables with eigenvectors |±, 0j, 6i) = cos(±| + Oi)\H) + sin(±| + 
6i)e^'^^\V), {i = 1,2) and eigenvalues ±1. The kets \V) and \H) represent here two 
orthogonal qubit states. One can write the annihilation operators Ai± of the photons in 
the modes observed by the detectors in terms of the annihilation operators in the input 
modes of the polarization analyzers. This gives 

Ai± = cos(±^ + ei)ai^H + sin(±^ + 9i)e''^'aiy 

A2± = cos(±^ + 92)a2,H + sin(±^ + 92)e"f''a2y. (3) 

We make the usual simplifying assumption that the probability P2{Diri , -D2r2 02, ^i, ^2) 
of detectors Di+ and D2+ to click is proportional to 

piD^r,,D2r,\<Pl,<P2,e,,e2) = {^{ua,^, U^.^J^) = {^\V^V\^) , (4) 

where V = A2T2^iri and ux is the photon number operator in mode X. This 
approximation is justified for linear detectors, and allows us to get simple analytic 
formulas estimating the detection frequencies. We obtain: 



P2 ~ p(L'lri,-D2r2l01,02,6'l,6'2) 

cosh^ K (I 1 



-1 + 3cosh2irl2 ^^^'^2 



sin(2^i)sin(2e: 



2 



}sinh^ K 
+ _i^3cosh2ir - 

For details of the calculation see Appendix A Note that for small values of K the 
probabilities approach the ones for an ideal maximally entangled two qubit state. 

As local detectors click at random, the formula for the "two-photon" interference 
visibility in this experiment can be put in a standard form 

max min -1 

V, ^ ^ = (6) 

where the maximum (minimum) is taken over (pa, (pt, da-, db- For small the visibility 
is close to 1 (see Fig. [2]). This case corresponds to a weak pumping field. In such a case 
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Figure 3. An experiment with many PDC sources (compare fig. 1). The symbols 
a.i (Ai) denote spatial modes before (after) polarizing beam splitters (PBSs). A^-i- are 
modes observed by detectors Di±, which are behind polarization analyzers PA^. 



a pump pulse usually produces nothing, and only sometimes a single pair. The prob- 
ability for other events is extremely low. If one increases the intensity of the pumping 
field, higher order (multifold) emissions start to play a role, and the visibility decreases. 
For the unphysical case of infinite K, the visibility is equal to 1/3. A classical radiation 
of a thermal type gives maximal visibility of this value in two-detector experiments. All 
this recovers the results of [TO] . 



2.3. Multi source case: GHZ-states 

Let us now extend the above scheme by introducing additional sources (see Fig. |3]). 
Consider an A^-photon experiment with N/2 sources. Modes n and (n + 1) meet at 
a polarizing beam splitter (PBS). Again, the detectors are placed behind two-channel 
polarization analyzers. A set-up of this kind was designed and used, for small pump 
power, to observe GHZ-type correlations [TTl 112] . 

The annihilation operators An± of the modes observed by detectors Dn± can be 
written in terms of the annihilation operators in the input modes of the polarization 
analyzers in the following way: 

^n± = ^(A^±e^'^Mr). (7) 

To simplify the discussion we assume that all 6'j's of the local settings are zero. In turn, 
modes A^^ and A]^ {n = 1, N) can be expressed in terms of the initial modes: 
aH _ H. aH _ 

- ■ - (81 
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Figure 4. Visibilities versus K for multisource GHZ-type experiments. The points on 
the plots denote critical values beyond which the visibility is not sufficient to violate 
standard Bell inequalities. 



where n = 1, — 1 and = for n = 1, A^. In these formulas we assume that 
H polarization is transmitted whereas V is reflected. 



Performing similar steps as in the previous case (see Appendix A ) we get the 
probability of a joint detection in the form of 

PN = p{Dir^, ...,DNrN\4>U ■■■,4>n) 



2 tanh^ K ^ JJ cos ^ 0^ j + 



Pn 



where 



N 



Pn = Y1 



j=0 



4J 



"i-i 



,fc=0 



tanh^^' K. 



(9) 



(10) 



The first term in the sum is equal to 1, while the last one reads 2^~Hanh^ K. In 
particular, 

P2 =1 + 2 tanh^ K, 

=1 + 8 tanh^ K + 8 tanh*^ K, 



1 + 18 tanh^ K + 48 tanh^ K + 32 tanh^ K 
ai 

+ 128tanh^K 



1 + 32 tanh^ K + 160 tanh^ K + 256 tanh^ K 



Pio = 1 + 50 tanh^ K + 400 tanh^ K + 1120 tanh^ K 
+ 1280 tanh^ K + 512 tanh^° K. 
The visibility reads 



N 



Pn 



Pn 



1 

Pn' 



:iii 



(12) 



Pn Pn 

In Fig. Hlwe compare the visibilities Vn as functions of K. For the unphysical very high 
values of K the visibilities asymptotically converge to the following value: 

2 



Vn{K 



oo 



(V2 + i)^ + (v^-i; 



N 



(13) 
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These values are related to modified Fibonacci numbers. Originally, they are 
elements of a sequence given by F{0) = 0, F{1) = 1 and F{n) = F{n — 1) + F{n — 2) 
thereafter. In this case, the sequence satisfies F'{0) = 0,F'(1) = l,F'{n) = 2F'{n — 
1) + -F'(n — 2). The solution reads 

F'(„) = (i + v^)"-(i-v^r, (14) 



Thus 



For example. 



V2 



y^iK ^ oo) = • (15) 



V2{K ^ oo) = 1/3 Vi{K ^ oo) = 1/17 

Vq{K ^ oo) = 1/99 Vs{K oo) = 1/577 

Viq{K ^ oo) = 1/3363. 

Note that for the analysis of these GHZ-like correlations we use only A^-photon 
interference described by the probability formula ([9]). Therefore, there is no ambiguity 
in the definition of Vat, as formula ([9]) does not reveal any interference involving less 
than photons. Simply, one can easily show that the marginals like 

do not depend on the phases. Had there been lower order interference processes, one 
could face difficulties in interpretation of the meaning of the parameter V . Thus far, 
these were fully resolved only in the two particle case [T3] . 

We are most interested in the critical values of Kent such that for K > Kcru, the 
visibility in the experiment is not sufficient to violate the standard Bell inequalities 
[m [HI [16]. These are given in Tab. [H We also give critical values of K, which are 
necessary to reveal entanglement in noisy GHZ states [TTf [18]. 

The values of Kcru subtly decrease with A^. This means that for this family of 
experiments, the critical pump amplitude is always at almost similar level. However, 
since the value of the K for which we observe high enough interference is bounded by 
Kcrit, the probability of an n = A^/2-pair emission, proportional to {^^^^^) , is also 
bounded from above. Thus, since one cannot arbitrarily increase the pump power, with 
increasing N one needs more and more time for the experiment. In Tab. [T] we also 
present the critical values for the probability pi that only a single pair is created by 
a given source, the probability that each source emits exactly one pair (this a kind of 
an optimal event for the studied interference experiment), probability 1 — p^ that one 
has at least one pair emitted by a given source, and the probability that each source 
emits at least one pair of photons. In realistic experiments, the visibility also depends 
on other factors (see, e.g., section [3]). Thus, parameter K must be much lower not to 
additionally reduce the observed non-classical interference. 

As shown above, statistics introduces a fundamental limit on the visibility in the 
multiphoton interferometric experiment, in which we use pairs of entangled photons 
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Table 1. Critical values of parameters corresponding to V'^"*, which is necessary to 
exclude local realistic (upper rows) or separable descriptions (lower rows). Respective 

symbols denote: 

1 — pn = tanh^ Kcrit ~ probability, that one, or more, pairs of photons are created at 
one source; 

Pi - probability of a single pair creation at one source; 

(1 — po)^/^ - probability, that each source emits at least one pair of photons; 
(pi)^/^ - probability, that each source emits exactly one pair of photons. 



]\T \/crit 1 n il r\^\^l'^ { -n 

V J^crit i- Pn Pi Pn) [Pi) 




2 ^ 0.4911 0.2071 0.1642 0.207106 0.164213 
i oo 1 1 


local realism 
entanglement 


4 ^ 0.4697 0.1918 0.1550 0.036790 0.024030 
1 1.0613 0.6180 0.2360 0.381985 0.055726 


local realism 
entanglement 


6 ^ 0.4404 0.1714 0.1420 0.005033 0.002864 
^ 0.9877 0.57218 0.2448 0.187258 0.014670 


local realism 
entanglement 


8 ^ 0.4232 0.1597 0.1342 0.000650 0.000324 

^ 0.9757 0.5613 0.2159 0.101100 0.003651 


local realism 

eiitanglenient 


10 0.4127 0.1527 0.1294 0.000083 0.000036 
0.9735 0.5629 0.2460 0.056493 0.000902 


local realism 
entanglement 


oo 0.3695 0.1250 0.10934 
0.9730 0.5625 0.2461 


local realism 
entanglement 



created in the parametric down conversion process. In order to obtain non-classical 
properties of down-converted photons the parameter Kcrit cannot be exceeded. It 
means that for experiments with many sources the probability of a joint emission, which 
supplies enough photons to observe a detection in each observation station, decreases 
exponentially with n and the experiment must take much more time. To put it short, 
one cannot expect too much progress with multiphoton interference techniques similar 
to the one described here, without a significant increase in the collection and detection 
efficiency. High intensity pumping is counterproductive. 

3. Distinguishability problems 

Interference is perfect only if one has perfect indistinguishability of various quantum 
processes that lead to the detection events. In an experiment a mode mismatch might 
appear for paths of propagation which are mixed (crossed) at a beamsplitter. This 
lowers the interference contrast. Such a case is described in the first subsection. The 
other cause of degradation of interference, in the case of multi-source PDC experiments, 
is the strong energy-time correlation within a down-converted pairs. Via this property 
one could in principle distinguish photons originating from different sources. The second 
subsection presents a discussion of a method of erasing such correlations with narrow 
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Figure 5. Entanglement Swapping Scheme. Source PDC(l) emits photons either to 
spatial modes ai and 61, or to Ai and Bi. PDC(2)-to mode pairs 02, &2 or A2,B2- 
The swapping occurs when detectors ii and 12 click simultaneously. Gray horizontal 
rectangles represent non-polarizing beam splitters, white diagonal are phase shifters, 
and black ones - perfect mirrors. 



spectral filters. 

3.1. The problem of mode matching 

In experiments like the ones considered above one has to overlap radiation from different 
sources. This, causes alignment problems, which may reduce the visibility even further. 
We present here the simplest case. We shall see that misalignment introduces an 
independent factor reducing the visibility. 

Consider the scheme presented in Figure |5l It consists of two type-I PDC sources 
which produce photons of a fixed polarization in pairs of spatial modes and bi, or 
Ai and Bi {i = 1,2). Mode 61 is crossed at a beam splitter with B2, and similarly 62 
is crossed with Bi. Behind these beam splitters we place two detectors, labeled ii and 
22- In the ideal case, if both of them click simultaneously, the two outer photons are 
entangled. We now want to investigate how the two-photon interference attributable 
to this entanglement swapping is decreased by the fact that behind the beam splitters 
photons from modes bi and Bi might be partially distinguishable from those from 62 
and B2, respectively. More precisely, a different signal arrives to a detector, depending 
whether the photon had been reflected or transmitted by the beam splitter. For 
example, the second source may produce photons with a slightly tilted polarization. 
Therefore, the annihilation operators related to the second crystal must be put as, e.g., 
X2 = cos 0X2 + sinaX^, {X = b,B). By we denote the radiation mode of source 
2 which is distinguishable behind the beamsplitter form the modes originating from the 
other source. 

Simple forms of the Hamiltonians describing the processes in the first and second 
source are given by: 

T-CpDCii) = ixi^ibi + AiBi) + h.c, 



Interference contrast in multi-source few photon optics 



10 



npDci2) = ^X(a2&2 + A^B'^) + h.c. (16) 
In the Heisenberg picture the operators evolve as follows: 

a2(T) = cosh fra2 — sinh fr(cos a 62 + sin a;(6^)^), 
A2{T) = coshKAa -sinhi^(cosa bI + sin a(5^)"^), 
b2{T) = cosh cos a (cos a 62 + sin a 62 ) 

— sinhi^ cosaflg! 

i?2(T) = cosh i^' cos a (cos a B2 + sin a B2) 

— sinhi^ cosaAg, 

6^(T) = cosh sin a (cos a 62 + sina 6^) 

— sinhi^sina a\, 

i?^(T) = cosh sin a (cos a B2 + sin a S^) 

— sinhi^'sina A2, (17) 



and 



X{T) = cosh KX - sinh KY\ (18) 

where X = ai,bi, Ai, Bi, and, respectively, F = bi,ai, Bi, Ai. X{T) denotes the 
annihilation operator in the X mode after the interaction time T. 

The probability of a detector click is proportional to the number of photons present 
in the mode monitored by detector, i.e. 

X (n,,(+)(T) +n,x(+)(r)) (n,,(+)(T) +n,x(+)(r)) :), (19) 

where n^iT) = x^(T)x(T) is the photon number operator in mode x after the interaction 
time T. Note that we sum over in principle distinguishable detection event, thus we have 
here + nf (recall that in our example x and X"*" represent orthogonal polarizations.) 

One can write the operators of the modes observed by the detectors in terms of the 
operators in the initial modes in the following way: 

A,(+) = -^(a,, + A,e*<^^) (x = l,2), 

^1 = ^(^1+ 62), i2 = ^(52 + &l), 

ti=b^, ii = Bl (20) 

(21) 

After some algebraic manipulations the final formula for the visibility of the 
detection rate f|T9|) is given by: 

Vi{K, a) = cos=^ a V^iK). (22) 

Figure M shows the values of K and a for which the state created in the entanglement 
swapping scheme can violate a CHSH inequality [19]. The border values of K correspond 
to the visibility equal to 1/ a/2. 
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Figure 6. The PDC efficiency parameters K and the ahgnment parameters a for 
which V > ^ in the entanglement swapping scheme. Only in the shaded region the 
CHSH inequahty can be violated. 



3.2. Partial distinguishahility of emission acts 

If one considers the multisource scheme presented in Figj3] one immediately sees that 
the tight frequency correlations of pairs of PDC photons originating from an individual 
act of emission may ruin our interference effects. The frequencies of the members of a 
pair (photons labeled here as 1 and 2) sum up with great accuracy to the frequency of 
the pump, ujp, 



This leads to strong correlations also in the time domain (such a phenomenon could 
be called "an Einstein- Podolsky- Rosen" effect"). Thus the pairs coming from different 
emission acts can be distinguishable in many ways - and interference suffers. This effect 
is independent of those studied in earlier Sections, where solely the influence of photon 
statistics was studied. It also requires a different calculational approach. As it was 
suggested in [2] a proper use of filters and detection gates leads to a high interference 
contrast. Later on it was shown in [22], [23] that the practical answer is to use pulsed 
pump and proper filters; this approach is used here. 

We shall use the following simplified description of a two-photon state emitted by 
one of the sources that makes through the filters: 



where ap{uJx) is an annihilation operator describing a photon mode of frequency uj and 
polarization P = H,V. The index x = 1,2 enumerates the directions, into which the 
photons are emitted. The symbol /(a;) denotes the spectral function of the filters, and 
g{ujQ) the spectral profile of the pump field. The operators satisfy the commutation 
relation [ap^Ux), a^p,{uj'^,)] = 6pp'6xx'S{i^ — ^')- Such a simplified description of radiation 



(23) 




(24) 
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(see e.g. Fearn and Loudon [2D]) works quite well in the case when the directions of 
emission are well defined. The function A represents the phase matching condition and 
will be replaced by 6{u!o — ui — UJ2) in our considerations. 
Next, if one puts the state into the form 

SUQ) = (25) 

(where Sl^ denotes the operator in formula (12^ ). one sees that its generalization to one 
emission per source of Fig. 3 is given by |\E') = S'|25']4...S'|r_^^|fi). This will be the 
starting point of our description. The probability that all N detectors register a count, 
respectively at times ti, t^r, is described by the usual formula for the iV-fold detection 
frequency: 

p(ti,...,t^) oc (M/|E(7)(ti)...i5i;)(t^)Ei+)(t^)...<)(ti)l^), (26) 
where the simplified effective field operators are given by 

= / rf^d,e-^-^.*'=a,,(.;,J, (27) 

with indices dk denoting both the beam leading to the detector and the polarization, as 
each detector observes only radiation of a specified polarization. The final annihilation 
operators a^^ are related to the initial ones by the usual algebra for the polarizing beam 
splitters and polarization analyzers. 

Following MoUow [2T], p{ti,...tN) can be written down using a square of the 
following amplitude: 

(0|Ei;Ht^)-<^(ti)l^)- (28) 
This is due to the fact that |\E') is an A^-photon state, therefore the action of N 
annihilation operators leaves out only the vacuum component. Such an amplitude will 
be denoted by ip(ti, t^). 

We assume that the detection stations measure elliptic polarizations, which can be 
represented by the following polarization mode transformation: 

aDiuj) = -^{(^hH + e'^avH) (29) 

(non polarization indices are dropped here). We assume that the polarizing beam 
splitters that mix the beams from two different sources transmit the H light and reflect 
the V light. 

An iV-photon interference is possible only if: either all PDC photons are 
transmitted, or all down-converted photons are reflected by the mixing polarizing beam 
splitters (this effectively leads to a GHZ state, {\H)^^ +\V)^^)/V2). The interference 
depends on indistinguishability of the two processes. Thus we have an overall amplitude: 

i^iti, ...,tN) oc ^pR{ti, ...,tN) + e'^ijTiti, tiv), (30) 
where e*^ represents the phases due to the measurement of the elliptic polarizations 
= ^t=i(pDi), and ipR and ipT are amplitudes with all photons reflected and 
transmitted, respectively. One has 

Pit,, ...,t^) cx IV^hI' + iV^Tp + 2Re [c'^^Pr^Pt] ■ (31) 
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It the above formula we have assumed that ■^'s are real, which will be the case for 
amplitudes used below. 

However, the detectors have a finite resolution time, which in comparison to the 
standard, in such experiments, pump pulse duration and the filter coherence time can 
be treated as infinitively long. Thus the overall detection probability behaves like 

/oo /'OO 
dti... dtNPiti,...,tN). (32) 
-oo J —oo 

Therefore, it is clear that the visibility is given by 



/oo POO 
dti... / dtNlpR{ti,t2,t3, ...,t]y_i,t]y) 
oo J —oo 



X 1pT{ti,t2,t3, ...,tN-l,tN) 



(33) 



Of course, above we have tacitly assumed normalization: 

/oo roo 
dti... dtN\Mtl,t2,t3,...,tN)\'' (34) 

-oo J —oo 

for X = R,T. 

It is important to notice that 

'ipniti, t2, tiv) = ^(tl, t2)...ij{tN-l,tN) (35) 
and, similarly, due to the action of PBS's, 

'0t(^1, ^2, ■-, ^JV) = V'(^l, h)'4^{h, ^5)^(^4, t7)...^{tN-2, tw), (36) 

where 

ijj{t,t') (X J dujQ j dui J duj2f{uji)f{uj2)g{ujQ) 

X A{u'q -uji- uj2)e'^^'e^^^' = 

= {m'Cilm^2{^)S\2\^) (37) 
(the polarization was taken here as H just for the sake of definiteness). 

If one assumes for simplicity that all relevant functions are Gaussian, that is 

and 

^(-) = exp I - V L (39) 



2ao 



one obtains that 



i){txM) = -^Q-^v{oix{t\ + t\) + a2txt2), (40) 



Interference contrast in multi-source few photon optics 



14 




Figure 7. The visibility of iV-particle GHZ-type interference produced in the setup 
of Fig. [3] The curves are drawn for = 4 (the highest), 6, 10 (the lowest). 



Note, that amplitudes of ijj, as promised, are real. 

We are now ready to calculate the visibility of A^-photon interference effects. For 
all photons reflected by the polarizing beamsplitters, the amplitude ipR takes the form 



I ( ^ 

i}R{ti,...,tN) = p]^exp ai^t- +a2^t2i-it 

\ i=l 1=1 

whereas, for all photons transmitted we have 

1 ^ 
'lpT{tl,...,tN) = fj^Gxp (ai'^t^i 



(43) 



1=1 



+ a2{tit3 + t2t5 + + - + tN~2tN)j ■ (44) 

Please note that the visibility is a multiple integral of a Gaussian function. The argument 
of the exponent in ipR{ti, ....ti^)%pT{tii ■■■■ti^) is at most quadratic in all t's. One can now 
extract the coefficients of a polynomial of ti and apply 



exp(—ax^ + l3x + 'y)dx = \ — exp I h 

V « V4a 



7 , (45) 



(a > 0). Calculations of such a kind have to be done times (see Appendix B). If one 
puts / = o-f/o-Q, the formulas for visibilities read: 

1 + f 



^4 = ^^^i^, (46) 



Va = (47) 

' (l + /V2)(l + 3/V2)' ^ ^ 

V 2(l + 2/r^ 

' (l + /2)(/2 + 2 + v^)(/2 + 2-v^)' 

Fio = 16(1 + 2f )V(5(/2 + 3 + v^)(f + 3 - v^) 

X (/2 + i + i/v/5)(f + l-l/v^)), (49) 
V,2 = 16(1 + 2f )5/V((l + f )(2 + f )(2 + 3f ) 

X (/2 + 2(2 + v^))(/2 + 2(2-v^))). (50) 
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Table 2. Values of /^'''* and for iV = 4, 6, 8, 10, 12, oo. 



TIT -rrtemp fcrit pcrit,approx 

^crit In In 



4 


1 

2^2 


3.806 


4 


6 


1 

A^/2 


3.592 


3.884 


8 


1 

8^2 


3.630 


4 


10 


1 

16^2 


3.700 


4.125 


12 


1 

32^2 


3.737 


4.237 


oo 





? 


4v^ 



For / <^ 1, that is narrow filters, V^*^™^ has the following tailor expansion: 

l^r^^l-^ + ^ + ..., (51) 



while for broad filters, / ^ 1, 

3iV^ 

Vn ^ ^^/"^^'- (52) 

The last equation allows us to easily find a approximate critical value of /, above 
which local realism certainly cannot be falsified (with the use of WWWZB inequalities 

Table [2] lists the critical visibility, required to violate the WWWZB inequalities, and 
the critical values of / at which we have this visibility obtained from (14611501) and flS^ . 

As can be seen from Table [21 jg upper bound for /^'*. The 

approximation becomes less and less accurate with growing A^. On the other hand, 
the sequence of /^** appears to be monotonously increasing for > 6. Hence we can 
conjecture that 3.737 < Z^'** < 4^2. 

4. Structure of Noise 

From ([5]) one can see that the effective, or "apparent" , two qubit state produced in the 
two-photon experiment of Fig. 1 is a Werner state (a mixture of a maximally entangled 
state 10+) = j^{\HH) + \VV)) and the "white" noise), 

Pe// = 1^10+) (0+1 + ^, (54) 

where V is given by ([6]) (we hope the usage of the same symbol for visibility and vertical 
polarization does not cause any trouble for a careful reader). By the effective state we 
shall mean an A^-qubit state endowed with set of probabilities proportional to the one 
estimated in the actual experiment, and which are afterward interpreted as "A^-photon" 
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(|)4.94 

Figure 8. The setup used in [TB] for generating correlations as of l^p^) state. PDC- 
nonlinear crystal, BS-beam splitter, PA-polarization analyzer. 

coincidences. Interestingly, situation like the one of eq. (IMI) is not necessarily the case in 
experiments with more observers; the noise can be, in general, structured (or "colored"). 

In this Section we aim to analyze the noise admixtures in the effective output state 
in the scheme of Ref. [16], presented in Fig. [HI The setup consists of a single source 
of entangled photons, two non-polarizing beam-splitters (BS), and four polarization 
analyzers. If the source produces just two entangled pairs, in the state -^{IHV) + \V H)) , 
and one registers one photon in each of the detector, this collapses the initial state onto 
the so called \ip4) state: 

IV'4) = ^/T/3{\HHVV) + \VVHH)) 

+ ^/T/12{\HVHV) + \HVVH) + \VHHV) + \VHVH)). (55) 

The first two positions in the kets correspond to side a in the figure, the other ones to 
side h. 

However, this is an idealized situation, in which we do not take into account higher 
order emission processes. As we shall see these processes may influence experimenter's 
estimate which state was actually produced in the experiment. These corrections, due 
to the statistical properties of the PDC radiation are unavoidable, as the efficiency of 
the PDC process must be quite high to have a significant probability of having of double 
pairs produced by a single pump pulse. This makes triple emissions also quite probable. 
Below we shall study the influence of these statistical effects on the interpretation of the 
experimental data as far as the estimate of the final state is concerned. Other effects 
may also influence such an estimate, but they will be ignored. 

The effective four qubit state observed in the experiment can be constructed via 
a tomographic method. From the probabilities of registering clicks in detectors we 
construct the correlation function, 

^(611, 01, ...,614,04) = ^ rir2r3r4 

f'l,»'2,''3,»'4=±l 

X P(^l,ri, -D2,r2, -Ds.ra, -D4,r4|01, 02, 03, 04, 6*1, 6*2, 6^3, 0^). (56) 
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The probabilities here are now normahzed, that is 

l,ri, -D2,r2, -Da.rg, -D4,r4|01, 02, 03, 04, ^1, 6*2, 6*3, ^4) = 1- (57) 

ri,r2,r3,r4=±l 

Note that the normahzation factor is the probabihty that one observes a chck at each 
of the four detection stations. 

Subsequently, the values of the correlation function for specific sets of angles give 
us the elements of the correlation tensor, 

01,fcl, ^2,fe2) 02,fc2> ^3,fc3> 03,fc3, ^4,fc4, 04,fc4), (58) 

where 

{0 for fci = 1 ( for ki = 1 

for A;i = 2 , 0,,^, = I 7i/2 for h = 2 (59) 
7r/4 for A;i = 3 [ for h = 3 

The three values of ki correspond, in the "spin picture" of a qubit, to three different 
complementary measurements associated with the Pauli matrices o'l = crx,cr2 = cry and 
= 0"^. We also need all marginal correlations, which are obtained from, for example, 

^(6*1, 01, 6^2, 02, 6^3, 03) = E rir2r3 

ri,r2,r3,r4=±l 

X p{Di^r^,D2,r2, D^^rs, -D4,r4l01, 02, 03, 04, Oi, 62, 63, 64), (60) 

That is we compute the marginals out of the observed four-fold coincidences. Please 
note that, the formula fl60p is not the one obtained in the experiment with observation 
three-fold coincidences in detection stations 1,2 and 3, even if the detection efficiency 
is perfect. All reasonings are based on the assumption of registration of a four fold 
coincidence (which collapses the initial state). The set of 256 numbers {Tijki}fj k^i=Q 
allows one to reproduce the effective density matrix: 

1 ^ 

P'^ff = 5Z '^ijki(Tl(T'^(Tlaf, (61) 

i,j,k,l=0 

where ctq = 1. This procedure leads to 

3 + 4cosh2i^ + 5cosh4/s: 
- 3(3-4cosh2ir + 5cosh4i^)'^^^^^^' 

1 -4cosh2ir + 5cosh4K A,,,,,, 



6(3-4 cosh 2K + 5 cosh AK) 

(62) 



i=l 

3-4cosh2fs: + cosh4/s: 
6(3-4 cosh 2K + 5 cosh AK) .^^ 



where are given by: 



#1 
#2 
#3 



^{\HVVV) + \VHVV) + \VVHV) + \VVVH)), 
^,{\HHVV) - \VVHH)), 

\{\HHHV) + \HHVH) + \HVHH) + \VHHH)), 
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Figure 9. The plot of the visibihty (sohd), the fidehty (dashed), and Vtotai (dotted) 
of Peff as functions of K. The vertical line represents the critical visibility for the 
violation of the standard Bell inequality. Parameter e is above 1 for all values of K. 



• #4: 


\VVVV), 




• #5: 


1{-\HVVV)- 


\VHVV) + \VVHV) + \VVVH)), 


• #6: 




\HVHV) - \HVVH) - \VHHV) - \VHVH) + \VVHH)) 


• #7: 


l(-\HHHV)- 


- \HHVH) + \HVHH) + \VHHH)), 


• #8: 


\HHHH). 





The usual definition of visibility, the one used in earlier sections for GHZ 
correlations, does not reflect the quality of interference in such an experiment see Fig 
[9l Therefore we present also the fidelity of pe// with respect to the desired \ip4) state, 

F = {ipilpeffllpA). 

A clearer picture of the interferometric properties of the state is given by the 
following function of purity 



total 



(63) 



where N is the number of detecting stations ("qubits"). Note that Vtotai varies between 
1 and 0. The square root takes into account that correlations tensor components Tij^i 
enter Trp^ squared. The parameters are presented in Fig. [9l Note that the formulas 
used for the visibility in the earlier section do not reflect here the real situation. It is 
clearly visible that Vtotai is a much better parameter reflecting the loss of interference in 
this case. 

The value of the parameter 

g ^ ^i,j,k,l ijkl .g^x 

maxii^ ' 

where, ma.xE is the maximal value of the correlation function fl56l) [18], for the critical 
value e = 1 gives one the threshold for having a sufficient condition for entanglement of 
the apparent state. Thus one could introduce yet another interferometric parameter of 
relevance. In the case studied here e > 1 for all values of K. That is, the apparent state 
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is always entangled. Finally, max^T^^-^^ with summations over pairs of orthogonal 
observable directions gives us the extent to which the given state may violate Bell 
inequalities [21]. These results will be presented elsewhere. 

Our discussion just a sketchy representation of the real situation. For real setups one 
must use a much more involved analysis, including a much more refined (and consistent) 
description of the detection process. Nevertheless, we think, the presented results signal 
specific problems that one might expect in real experiments. 

5. Conclusions 

We have reviewed some difficulties that might be encountered in a typical few-photon 
PDC experiment. The possibility of creating more than the required number of pair of 
entangled photons in the parametric down- conversion, or having some leftover frequency 
correlation between photons of the same pair, is unavoidable from the very fundamental 
point of view. 

There are some fundamental limitations for a successful multiphoton interference 
PDC experiment. The maximal value of the process efficiency parameter allowing to 
violate a standard Bell inequality decreases from K{N = 4) = 0.4911 to K{N oo) = 
0.3695. This makes the production of at least one pair in each source exponentially 
improbable with a growing A^. Hence, as a function of A^, the whole experiment requires 
at least exponentially many runs, not only due to the growing number of measurements. 

As pointed out in [21 [221 [23], the frequency correlations the photons can be removed 
by better defining their frequencies with narrow spectral filters. By the Heisenberg 
uncertainty principle, this causes the instants, at which the photons of a pair reach 
the detectors, to be sufficiently undefined to wash out the complementary temporal 
correlations. We have shown that one for pulse pumped sources cannot violate a 
standard Bell inequality for crj/ao > A\/2. 

We have also shown that in the case of some experiments, the noise introduced by 
unwanted additional emissions of photon pairs may have a quite complicated structure, 
which may affect the interpretation of the experimental results. 

The description that we used was as simple as possible. We conjecture that 
refinements would lead to even more pronounced effects that those shown here. For 
specific experimental setups one must perform a more detailed analysis, involving a 
more realistic description of detectors response, and their inefficiency. 
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Appendix A. Derivation of formula ( 15]) 

This appendix shows the way to calculate the formula ([5]) for p{Dirj^, D2r2\4>i, 02, Oi, Oi)- 
Using ([2]) and ([3]) we obtain 

^ oo 

n,m=U 

( H 71 

X (^cos(6'i + ri-) cos(^2 + ?"2-)^ 

X \n - l)^,\m)l\n - l)l\m)l 

+ cos(6'i + ri—) sin(6'2 + r2—)y/nm e^'^' 

X \n - l)^,\m)l\n)Z\m - 1)1 

+ sin(6'i + T'l—) cos(^^2 + r2—)\/nm e*"^^ 

X - l)l\n - l)Z\m)l 

+ sin(^i + ri^) sin(^2 + r2^)m e'^'^'^'*'^^ 

X i<i"^-i)rj<i^-i)r.) (A.i) 

and the probability 

D2r, 101, 02, 9u 62) = U-^\Vm'' 

^ 00 00 

= Af-' y y tanh"+'"+"'+"^' 

cosh^K ,Vn 

n,m=0 n' ,m'=Q 
^ (^^n—l,n' — l^m,m'^n—l,n' — l^m,m' 

X cos^(6'i + ri—) cos^(6'2 + ^2^) nn' 

X cos^(^i + '"i^) sin^(6'2 + r2 — )Vnmn'm' 

~l~ '^n,n"^m— l,m'--l'^n— l,n' — l'^m,m' 

-)cos (^2 + r2- 



X sin^ (^^1 + r 1 — ) cos^ (6^2 + r2 — ) V nmn'm' 



IT TT 

X sin^(0i + '"i — ) sin^(02 + r2—)mm' 

~l~ '^^n—l,n'^m,m' — l^n—l,n'^m,m' — l 

X sin(2^i + ri|) sin(2e2 + r2^)W e^(*^+'^^) 

~l~ '^^n' —l,n^m' ,m—l^n' — l,n^m' ,m—l 
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X sin(2^i + ri|) sin(2^2 + r2^)nm e''^'^'^'^^'^ 



(A.2) 



After the summation over n' and m' many terms vanish and we get the final formula 
for the probability: 

p{Dir^, D2r2\(pl, 4>2, Oi, 62) 



P2 



cosh* K 



J2 tanh2("+'") K 



n,m=0 



X 



I cos^(ei + ri^) cos2(e2 + r2^] 



+ sin 



cos 



h + ri-j sm [62 + r2-) m 



h + ri-j sm [62 + r2-) nm 
+ sin^(6'i + '"i^) cos^(6'2 + nm 

1 TT TT 

+ - sin(26'i + ri-) sin(26'2 + r2-) n{m 



IT 

COS^(6'i + '"l^) COS^(^2 + ^2 



1) COS I 



4' 



+ sin 



ri — sm 

A' 



h + r2 



A' 



+ i sin(2^i + n^) sin(2^2 + ra^) cos (0i + ^2) 



+ tanh^'K 

,2 



+ sm 



+ cos 



cos 



h + ri-) cos 



TT, 



r2- 



'1 + ri-j sm" 

1 ^\ • 2 
^1 + ri-) sm 



?2 +^^2 



5*2 + ?^2 



+ sin2(^i +ri-)cos^(^2 



-) 

vr, 
4' 



?^2 



4 'J 



(A.3) 



Note that the terms in the bracket after "tanh^ fT" sum up to 1 and Af = 
(-i+3cosh2X)^tanh K rjj^^ valuc of J\f wc get usiug the normalization equation, namely 

Sri r2=±i^'(-^ifi' -^2r2l0i5 02, 6*2) ~ ^1113'%) usiug the above facts, and standard 
trigonometric relations, one can rewrite Eq. (IA.3P to get the required final formula. 
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By eq. f l33l) . the integral of our interest reads: 



Va 



f-2 I dt. 



dt-2 



oo 



oo 



dtA 



X exp 



2aj 



3 / "''4 

oo J — oo 
4 



^2 + 2a' 



The argument of the exponent can be put as 



(B.l) 



^a2 + 2ar' 



7 



f 
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(B.2) 



We have singled out the part of the expression which depends on t\. After integration 
over ti, which employs (l45l) . we obtain 



7r 



r2j2a2 



dtn 



dt. 



did 



X exp 



2a} 



i=2 



(^'jt4:(t2 + ts 



{t2 + t3f(yj 

al + a) 



(B.3) 



Again, we define the following quadratic form /2(t4), in which t^ is singled out: 



2a 



al + a^ 



f ^2 



2ap 



ti + 2 



{t2 + hVj 

al + a) 



a 



^-2{t, + t,fa} + {tl + tl){al + a)f 
f 



2{a^ + Sa^a] + 2a j) 
and integrate in fIB.SP over td. This brings our integral to the form 



(B.4) 



P(2a^- 



dto 



dt-i exp 



K + af)K + 2af) 



2a|Kfe + t3)^ + K + af)^(ti + ti 
a^ + 3aoVJ + 2aJ 
We put the function in the exponent into the following form 



(B.5) 



2{al + 2alal + 3a 



7 



af, + ?,ala) + 2a^ 



^a^t^ 
aj-t^ 



r^3 



+ Sa^aj + 2aj 
2{a', + 2alal + 3aj) 
a^ + 3aiaj + 2aj 2- 



(B.6) 



Therefore, the integral reduces to: 



2ajr 



y ((t2 + aj){a^ + 2ala} + 3a^) 



dt2 exp — 2(t1 



^0 + 



/ +2 



^a2 + 2af2 
and after purely algebraic simplifications gives 



(B.7) 
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